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Abstract 

We introduce a novel method to investigate the stability of wave packet dynam- 
ics under perturbations of the Hamiltonian. Our approach relies on semiclassical 
approximations, but is non-perturbative. Two separate contributions to the quan- 
tum fidelity are identified: one factor derives from the dispersion of the wave pack- 
ets, whereas the other factor is determined by the separation of a trajectory of the 
perturbed classical system away from a corresponding unperturbed trajectory. We 
furthermore estimate both contributions in terms of classical Lyapunov exponents 
and find a decay of fidelity that is, generically, at least exponential, but may also be 
doubly exponential. The latter case is shown to be realised for inverted harmonic 
oscillators. 
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1 Introduction 



It has long been appreciated that, in contrast to chaotic classical dynamics, the time 
evolution of a quantum system shows no sensible dependence on initial conditions. This 
follows immediately from the unitarity of quantum dynamics. Hence, rather than being 
concerned with instabilities at large times, the notion of quantum chaos is commonly 
reserved for semiclassical studies that aim at relating statistical properties of stationary 
quantum states to dynamical properties of chaotic classical systems (see, e.g., [HE])- More 
recently, however, the behaviour of quantum time evolutions at large times has attracted an 
increasing attention (see, e.g., 0111 EJE])- Both in the dynamics of observables (Heisenberg 
picture) and in the evolution of wave functions (Schrodinger picture) it has been proven 
that there exists a time scale (the so-called Ehrenfest time) , depending on the semiclassical 
parameter h, below which the quantum dynamics can be well approximated in terms of 
the associated classical time evolution. Moreover, if the classical dynamics are chaotic, this 
time scale is inversely proportional to a suitable classical Lyapunov exponent. 

Some time ago Peres suggested [7j that instead of studying the behaviour of quantum 
dynamics under a change of initial conditions one should investigate its stability under 
perturbations of the Hamiltonian: Suppose that an initial state ip is evolved under the 
unitary dynamics Uo(t) generated by the quantum Hamiltonian Hq, one compares this 
with the evolution U £ (t)ip determined by the perturbed Hamiltonian H £ = Hq + eV . Here 
V is a perturbation of unit strength, and e is a variable parameter. Then the quantum 



measures how sensibly the dynamics reacts to this perturbation. It can also be viewed 
as a means to quantify to what extent the initial state can be recovered after it has been 
propagated for a time t with the unperturbed dynamics, and then the time evolution is 
reversed with a perturbation being turned on. For that reason the quantity (jl.lj) is also 
known as quantum Loschmidt echo. 

Peres analysed F(t) in perturbation theory, and found an initial decay F(t) ~ 1 — 
Cy ^(et/h) 2 . Since the reliability of perturbative results requires et/h to be small, one could 
view Peres' result as indicating a Gaussian decay of the fidelity on an initial time scale that 
depends on e and h. Later work focused on time scales beyond this perturbative regime or 
on strong perturbations, respectively, and found an exponential decay |H1 HE EH UJ] • Its rate 
is either determined by Fermi's golden rule {8, or, for stronger perturbations, by a classical 
Lyapunov exponent [HUH]. Further studies related the behaviour of the fidelity to the decay 
of (quantum as well as classical) correlations ^2] ■ All of these results rest on a number of 
approximations and assumptions. Hence the precise time scales for the different regimes 
depend on various factors as, e.g., initial states, strength of perturbation, averages over 
random perturbations, and dynamical properties of the corresponding classical dynamics. 

Here our principal aim is to develop an alternative approach to the decay of quantum 
fidelity for particular initial states. The method that we shall introduce below is non- 
perturbative (quantum mechanically as well as classically) and employs only semiclassical 



fidelity 



F(t) = \^A{tr i mm 2 




2 



approximations with a rigorous control over the errors. Previous semiclassical studies of 
fidelity decay used the Van Vleck-Gutzwiller propagator for the time evolution of Gaussian 
initial states [01 E3 El 03] • This procedure takes care of the leading term in an expansion 
in powers of H, with an error that is, formally, smaller by a factor of h. For finite times this 
is indeed true, but there is no analytical control over the semiclassical error that arises in 
estimates of the fidelity decay at large times. Our method, however, allows to bound the 
semiclassical error in terms of the linear stability of an associated classical dynamics. It 
can in particular be applied to Gaussian states. And although our approach requires no 
particular assumptions about the nature of the classical dynamics, we are mostly interested 
in the case of chaotic (i.e., exponentially unstable) classical trajectories. In that case we 
find a decay of fidelity prior to the Ehrenfest time that generically is at least exponential. 
But it may also be doubly exponential, which we show to be the case in the example of 
inverted harmonic oscillators. 

This paper is organised as follows: In Section El we introduce the wave packets that 
we shall consider as initial states and review their semiclassical dynamics. The decay of 
quantum fidelity is investigated in Sectional with an emphasis on the behaviour of Gaussian 
states. An exact calculation of the fidelity for inverted harmonic oscillators is performed in 
Section |3J We finally summarise our findings in Section El Three appendices are devoted 
to a number of technical considerations: We first review the metaplectic representation, 
then discuss a transformation of positive-definite, symmetric matrices to a diagonal form, 
and finally collect estimates of matrix norms and singular values. 



2 Localised wave packets 

The initial states to which our approach applies are wave packets with a localisation both 
in position and momentum. By this we understand a concentration of the quantum state 
in a suitable phase space representation on a single point, when the semiclassical limit 
is performed by passing to a small (effective) Planck's constant h. We specify the wave 
packets in terms of normalised, smooth, and rapidly decreasing functions 4>(x) (Schwartz 
test functions) of x e M. d . Examples for this are provided by the Gaussian functions 



^(x) = ( det l™ Z ) 1/ \^ z *, (2.1) 



where Z is a complex, symmetric d x d matrix with positive-definite imaginary part. 
For the purpose of semiclassical investigations we introduce the scaling 

<p\x) = H- d ^<f>(x/Vn) . (2.2) 

This produces quantum states that are semiclassically concentrated at zero in position and 
in momentum. Such a phase space localisation is best analysed in the Wigner representa- 
tion 

W[<f) h }{& x) = J (j) h {x - f) <f) h (x + f ) e~i v * dy (2.3) 
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which, if multiplied by (2iih) d , converges to 5(£, x) as h — > 0. A subsequent application 
of the phase space translation 

D{p, q) = e -i(* p -P-Q) (2.4) 

therefore yields a wave packet 

4>f Ptq) (x) = e-^D(p, q )^(x) 
= e i p - {x - q) 4> h (x - q) 

with phase space localisation at the point (p,q). The phase convention made here is 
introduced for convenience; it merely simplifies some of the expressions below. 

The time evolution of such a state, generated by a quantum Hamiltonian H that arises 
as a Weyl quantisation of a classical Hamiltonian H(p, q), 

mx) =JjH(p,s?) e*^-») <%) ^ , (2.6) 

can be determined semiclassically [TJJ |3] to be 

U(t)<Pl, q) = e^Dfa, q t ) KSt) <f> H + O t (VK) . (2.7) 

The main term on the r.h.s. is again a wave packet of the type ([2.5)1 . but now localised at 
(p t , q t ). This is the point on the trajectory emerging in time t from the initial point (p, q) 
under the classical dynamics generated by the Hamiltonian H(p,q). Moreover, 

Rt= [ {q s -p s -H(p s ,q s )) ds (2.8) 
Jo 

is the action of this trajectory and St is the associated stability matrix. This is a real, 
symplectic 2d x 2d matrix that arises as a solution of 

S t = JH"(p t ,q t )S t , S = t. (2.9) 

Here J is the symplectic unit (|A.3|) and H"(p, q) is the Hessian matrix of the Hamiltonian. 
Equivalently, the stability matrix is given as 

(dPt dPt\ 
al al) ■ (2.10) 
dp dq J 

The wave packet at time t on the r.h.s. of ()2.7|) arises from the initial state <ft h through 
the application of a unitary operator consisting of two contributions: the first factor is the 
metaplectic operator fi(S t ) that provides a double valued representation of the symplectic 
group (of linear canonical transformations), see Appendix lAl and |T5 ^ I16 [ IT7j for details. 
As can be drawn from (|A.8|) . a metaplectic operator does not change the semiclassical 
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(2.5) 



phase space localisation. In (|2.7j) it is rather responsible for the dispersion of the wave 
packet. The second factor, D{p t , q t ), then provides a translation of the wave packet along 
the classical trajectory. Finally, the error term Ot(Vh) stands for a vector whose norm can 
be estimated from above by K(t) y/h. The function K(t) > contains the linear stability 
of the trajectory (p t ,q t ). If the latter is exponentially unstable with maximal Lyapunov 
exponent A > 0, the function K(t) grows like te 3Xt as t — > oo Therefore, as long as 
t <C Te(K), with an Ehrenfest time Te{K) = | logfi|/6A, the error term remains small. In 
contrast, if the trajectory is stable (in an integrable system or on a KAM-torus) the growth 
of K{t) is algebraic (like t 4 ) and hence Te(K) = C fir 1 / 8 . In any case, this finding enables 
one to extend the validity of the semi classical evolution (|2.7j) to infinite times, when h — > 0. 
We remark that the main term in (|2.7j) actually is the leading contribution in a systematic 
semiclassical expansion j^J. If this is carried on to the iV-th term, the error is O t (h N ^ 2 ) and 
can also be controlled up to T E {h). 



3 Fidelity decay 

The quantum fidelity of an initial wave packet of the type ()2.5j) can most conveniently be 
calculated in the Wigner representation, 



F(t) = \mWU)MW( P , q) )\ 2 =11 w\u (t)4^]{M W[(j £ {t)^l, q) \{M §- 



dx 



(3.1) 

We now introduce the semiclassical result (|2.7|l for the perturbed and for the unperturbed 
time evolution, respectively, to this expression. The corresponding unperturbed and per- 
turbed classical trajectories shall be denoted as q®) and (p%,qf). Exploiting the be- 
haviour 

W[D(p, q)4](t, x) = W[1>](t ~P,x-q) (3.2) 

of a quantum state ip in the Wigner representation under phase space translations, one may 
change variables and define (Sp t ,Sq t ) = (p° —pf,q^ — qf). For the leading semiclassical 
contribution one thus obtains 

F scl (t) = Jfw\ii(^^-8p t ,x-6q^W\fi(SI)^(M 7^ • ( 3 - 3 ) 

Since semiclassically the Wigner representations of localised wave packets approach 8- 
functions, after a division by (2nh) d the result ()3.3j) can be viewed as a smeared out 
classical fidelity that measures the separation (6p t , 5q t ) of the perturbed trajectory from 
the unperturbed one. 

For a more detailed study of the expression ()3.3|) we now restrict ourselves to Gaussian 
initial states of the form (j2.1j) with the scaling (|2.2j) . The action of a metaplectic operator 
on such states can be calculated explicitly [TH | fTH ^ fTT] . 

jl(S)(f) Z ' h = e^ a 4) s[z]fi . (3.4) 
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Here S[Z] denotes a map, given by the symplectic matrix S, on the space of complex, 
symmetric matrices with positive-definite imaginary part to itself, 



S[Z] = (AZ + B)(CZ + D)- 1 , S=(fj . (3.5) 



Furhermore, a is a Maslov phase defined through 

detlmZ \i/4 n N\-i/2 



The Wigner transform of such a Gaussian state is well known to be a Gaussian in phase 
space, 

W[(j) z ' n }{tx) = 2 <i e-*«' a 0' G *«'"0 , (3.7) 

where 

/ (ImZ)- 1 -(ImZy^ReZ \ 

z \-ReZilmZ)- 1 Im Z + Re Z(lm Z)~ x Re Z J 1 ; 

is a symmetric, symplectic, and positive-definite 2d x 2d matrix with unit determinant. 
The behaviour of ()3.8)1 under the transformation ()3.5j) can be inferred from an application 
of a metaplectic operator to a Gaussian state in the Wigner representation ()3.7|) . This way, 
from ()3.4|) and (jA.8|) one concludes that 

G s[ z] = (S- 1 ) T G z S- 1 . (3.9) 

Now, ([3.3)) is a Gaussian integral that can immediately be evaluated, and the result 
may be factorised according to 

F scl {t) = F disp (t) F class (t) . (3.10) 

The first factor 

-1/2 



F disp (t) = (det(G s o [z] + G sm )) (3.11 



is determined by the dispersion of the wave packets. This interpretation follows from the 
fact that setting 5p t and 5q t to zero in ()3.3)1 . and therefore removing the phase space 
translations that arise from ([2.7)1 . the integral would exactly yield ([3.1 lj) . In fact, F disp (t) 
measures the differences in the dispersions caused by the two dynamics in question. This 
contribution is independent of H. The time dependence of ([3.11)1 follows from the relation 
(13. 9 j) with S® and Sf, respectively. It is therefore completely determined by the linear 
stabilities of the perturbed and the unperturbed classical trajectory. In addition to this, the 
second factor F dass (t) is influenced by the actual separation (Sp t , Sq t ) of these trajectories. 
It reads 

F class (t) = 2 d exp \-\{5p u Sq t ) ■ G s?[z] (t - T^) (Sp t , 5q t )] , T t>£ = 1 + G^ [Z] G S![Z] 



(3.12) 
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and, despite of its ft- dependence, essentially represents a classical fidelity since it is localised 
on the separation of the classical trajectories: If one divides F(t) by (2nh) d as discussed 
above, the contributions of (j3.11|) and (j3.12j) converge to 6(5p t , 5q t ) as h — > 0. This also 
explains the necessity of H in ()3.12j) . We remark that expressions equivalent to ()3.1())) - (j3.12j) 
have independently been obtained by M. Combescure and D. Robert [18J. 

At this point we stress that in general the separation (Sp t ,5q t ) of the trajectories for 
large t differs essentially from the corresponding behaviour of the linearised dynamics. In 
particular, an exponential instability expressed in terms of positive Lyapunov exponents 
does not imply an exponential growth of the norm of (Sp t , Sq t ). In fact, for the dynamics 
of a bound system this quantity obviously is bounded. But even then the exponent in 
()3.12|) will often grow exponentially due to the presence of the stability matrices S® and 

si 

The contributions of F disp (t) and of -FciassW to the decay of fidelity will now be studied 
separately. This procedure makes sense if h simultaneously approaches zero in order to 
maintain the condition t <C Te{K). In this regime the individual contributions to F(t) 
determine the leading behaviour of the fidelity as t — > oo and h — > completely. 



3.1 Contribution of wave packet dispersion 

We begin with discussing the behaviour of F^ Bp (t) as t — > oo. Since both G s o^ and Gs%[z] 
are symmetric and positive-definite, we can convert these matrices into a diagonal form 
as explained in Appendix [Bl This implies that there exists a real, invertible matrix M t 
such that MfG S 0{ Z ]M t = 1, and at the same time MfGsnz]M t = D t is diagonal, with the 
(positive) eigenvalues A k (t) of 

G s°[z] Gs U z ] = S t G z 1 (( S ty 1 St) TG z(Siy 1 (3.13) 

on the diagonal. Furthermore, since Gz is symmetric and positive-definite, it is a square 
of a symmetric and positive-definite matrix, Gz = 7 2 . Hence, (j3.13|) is conjugate to a 
matrix NfN t , with N t = 7 (S^) -1 ^ 7 -1 . This means that the eigenvalues A^(t) of ()3.13|) 
are squares of the singular values of N t (see Appendix IH} • Since 7 is independent 

of t the time dependence of fik{t) therefore is asymptotically determined by the singular 
values p>k(t) of (5'f)~ 1 S' t . More precisely, there exist constants Ci,C 2 > such that 

CiMO < M*) < c 2 h(t) ■ (3.14) 

We also exploit the fact that the product (5'f) _1 5'° of two symplectic matrices is again 
symplectic. This implies that its singular values /ifc(t) arise in pairs of mutually inverse 
numbers. Thus, they can be ordered as in (jC.7|) . 

The determinant that yields -Fdisp(^) according to ()3.11|) can be evaluated as in Ap- 
pendix El see (IB. 1)1 . Taking into account that G s o^ has unit determinant and that the 
eigenvalues Ak(t) of ()3.13|) are given by squares of the singular values yUfc(t), we obtain 

-1/2 



-Fdisp(t) 



2d 



n( i+ M*) 2 ) 



Jk=l 



(3.15) 
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The estimates (J3.14J) then allow to bound ()3.15j) from below and above in terms of 



d 



-1/2 



(3.16) 



Y[^ k (t) 2 + 2 + j2 k (t)- 2 ) 

.k=l 

More specifically, there exist constants C 3 , C 4 > such that 

d d 

C 3 ll^ity 1 < Fdisp(t) < C 4 HM)- 1 . (3.17) 

k=l k=l 

Since the product over the inverse singular values contains only factors with fik{t) > 1, one 
can introduce the simple estimate 

d 

^max(t) < Y[fa{t) < /imax(t)^ 1 W) ■ ( 3 ' 18 ) 
k=l 

The quantities ftk{t) are singular values of a product of two symplectic matrices to which 
the inequalities (jC8|) may be applied. Thus, when choosing k — 1 in (jC.8|) . the l.h.s. of 
()3.18|) can be bounded from below by 

max ^ ToTT ' TFnT f < A^max^) • (3.19) 



/^max(»S'f ) A i max(*S't ) ) 

Furthermore, the maximal singular values of S% and 5*° determine the maximal Lyapunov 
exponents according to 

\°/ £ = limsup - log H^^Hhs 

J (3.20) 
= limsup - log/Vax^t ) , 

so that in case 5\ = A £ — A ^ the l.h.s. of ()3.19j) is asymptotic to e'^'* as t — > oo. 

The r.h.s. of 1)3.18)1 may now be estimated in a similar manner: Apply the rightmost 
inequality in (jC.8j) to each factor, and for the term with k = d use that fid{S^ £ ) = 1, see 
Appendix O This finally yields the upper bound 

[/imax(5'°) /imax(SI)] d 1 ' min{/i max (S'°) , /i ma x(SI) } (3.21) 

for /i ma x(^) d_1 P'dit) ■ Asymptotically, as t — > oo this approaches exp{[(d — 1)(A° + A £ ) + 
min{A°,A £ }]t}. 

The above estimates can be summarised to provide the following statement about the 
asymptotic behaviour of -Fdi sp (t): There exist constants C^,Cq > such that 

C s e- Lt <F 6iep (t)<C 6 e- Lt , (3.22) 

with 

\5X\ < L < (d - 1)(A° + A £ ) + min{A°, A £ } . (3.23) 

Thus, once the maximal Lyapunov exponent of the perturbed dynamics differs from the 
unperturbed one, the asymptotic decay of F disp (t) is essentially exponential. 
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3.2 Contribution of classical trajectories 

The remaining factor -F c i ass (t) that determines the decay of fidelity is influenced by the 
linear stabilities of the perturbed and of the unperturbed classical trajectories as well as 
by the separation (Sp t , 5q t ) of the trajectories. The contribution of the stabilities can be 
treated in a similar manner as above, whereas the behaviour of the separation is largely 
unknown in a general linearly unstable system. Precise estimates are rare, but can possibly 
be achieved in particular cases (see, e.g., Section |U and [T§]). 

A first simplification of the expression (j3.12)l can be achieved by introducing the matrix 
M t that converts Ggo[z\ and Gs*\z\ into a diagonal form. The exponent of (|3.12j) . without 
the factor —1/h, then reads 

M t -\5p t ,5q t ) ■ (1 - (1 + D t )~ 1 )Mf 1 (5p t , 5q t ) (3.24) 

where, as above, D t = G s e [Z]M t is the diagonal matrix with the eigenvalues A&(t) > 
on its diagonal. The quadratic form 1 — (1 + Dt)' 1 defined by (|3.24jl is positive-definite; 
its eigenvalues are A k (t)/(1 + A k (t)) > 0. Thus, whatever value the separation (5p t ,5q t ) 
attains, one immediately concludes that -F c i ass (t) < 2 d . And although the eigenvalues A k (t) 
may eventually grow as t — > oo, this quadratic form remains bounded. Any influence of the 
linear stabilities on -F c i ass (t) hence is encoded in the matrices M t as they appear in (|3.24jl . 

An upper bound for the expression ()3.24j) follows from the estimate (jC9|) derived in 
Appendix IU1 when choosing A = 1 — (1 + D t )~ 1 and B = Mf 1 . For the calculation of ||A|| tr 
we notice that this is given by the sum of the eigenvalues Afc(£)/(1 + A^(t)). These can 
be grouped in pairs with Afc(t) and Afc(t) -1 since the latter are eigenvalues of a symplectic 
matrix. Thus 

Afc(t) Akjty 1 

1 + A k (t) 1 + A fc (t)- 



d 



\A\L = > —pyr + K l > r - = d . (3.25) 

1 l|tr ^ UA, W UA, / -I 1 v ' 



k=l 

Furthermore, we observe that 

llMr 1 !!^ = HMtMjy 1 = \\G s?[z] \\ tT . (3.26) 
Using (j3.9|) . the rightmost expression can be factorised with the help of (j(].4|l . leading to 

l|Mr 1 ||L<ll^°llHsl|G'z||HS. (3.27) 
Our final upper bound for ()3.24j) therefore reads 

rf||G'z||Hs||^ ||^ s (^ 2 + ^ 2 ) . (3.28) 
For the contribution (j3.12j) to the fidelity this provides us with a lower bound of the form 

F clas& {t) > 2 d exp[-^||G z || H s ll^llls (*P 2 t + ^)] • ( 3 - 29 ) 

In addition, if the unperturbed trajectory possesses a positive maximal Lyapunov exponent, 
the factor HS^Hhs g rows asymptotically like e 2A °*, see ()3.20|) . 
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In order to achieve a lower bound for (|3.24j) according to (jC.10|) we first notice that the 
symplecticity of (|3.13|) implies A min (£) = A max (t) _1 = /i max (i)~ 2 . This leads to 

Amin(A) = — 1 > ~ /i max (t)- 2 . (3.30) 

1 T ^ 



Then ()3.14j) and (jC.6j) yield the further bound 



A min (A) > K^^iS^S?)- 2 

> K\ // naax (5'f ) A i max(*S't ) ) 

with some constant K\ > 0. We now require a lower bound for fi2d(B) 2 = /i 2 d(M t -1 ) 2 , and 
first notice that this quantity is the lowest eigenvalue of (MtMf)' 1 = G s o^, which in turn 
is the inverse of the largest eigenvalue of this matrix. Making use of the relation ()C.6j) we 
then conclude that 



A max (G ! 5 t [2]) = Vmax(((St) l ) T Gz{St) *) 
_ /'max (Gz) /'max 



(3.32) 



Collecting the above estimates therefore provides us with the lower bound 

^2/imax(^)"Vmax(^°)' 4 ((5p 2 + 5q 2 ) (3.33) 

for ()3.24|) . with some K2 > 0. Hence, 

^dass(t) < 2 exp — — — . (3.34) 

L pmaxi'-'t J h L ma,x\ J t J J 

Furthermore, in the case of linearly unstable trajectories (I3.20|) implies that 

/w(^) 2 /w(^ ) 4 ~ exp[(2A £ + 4A°)t] , (3.35) 

as t — > 00. 

We have so far refrained from estimating the squared distance 5p 2 + 5q 2 the perturbed 
trajectory can separate itself away from the unperturbed one. However, any further state- 
ments about the decay of -F c i ass (t) require some knowledge of the behaviour of that distance. 
But this, as already mentioned, seems to be difficult to obtain. E.g., one can in general 
not exclude that this quantity vanishes infinitely often (see [HI] for a particular case), or 
asymptotically approaches zero. Such situations may occur, if the perturbation V(p, q) 
is confined to a bounded part of phase space, but the trajectories are forced to leave this 
domain in a particular channel. According to ()3.12j) . at those instances where 5pf + 5q 2 
vanishes, -F c iass(£) clearly acquires its maximal possible value 2 d . On the other hand, if the 
energy shells corresponding to both the perturbed and the unperturbed classical Hamilto- 
nians are bounded, 5p 2 + 5q 2 will necessarily be bounded, too. 
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However, if the classical motion is unbounded (as in the example in Section^}, an upper 
bound for this distance would be helpful. In order to achieve such an estimate we consider 
a Taylor expansion of qf) about e = with remainder term of first order, i.e., 



(3.36) 



e'=9e 



where 9 G [0,1]. The derivative on the r.h.s. can now be identified as a solution of a 
differential equation in the variable t. Abbreviating 

z E (t) = ^- £ (P e t,Q e t), (3-37) 

the fact that (pq, q$) = (p, q) for all e implies the initial condition z E (0) = 0. Moreover, a 
derivative of Hamilton's equations of motion 

(plql) = JH' £ (pl,ql) (3.38) 

w.r.t. £ yields 

z £ {t) = JH':(plq\) z £ (t) + JV'(plq £ t ) , (3.39) 

where V'(p, q) denotes the gradient of the function V(p, q) whose Weyl quantisation yields 
the perturbation V of the quantum Hamiltonian. A solution of the inhomogeneous differ- 
ential equation (j3.39J) with the prescribed initial condition is then provided by the integral 

z £ {t) = St f(S?)- x JV'{jf„ q s s ) ds . (3.40) 
J o 

Used on the r.h.s. of ()3.36|) this expression allows us to relate the separation (5p t , 8q t ) of 
the trajectories to their linear stabilities and properties of the derivative of the classical 
perturbation V. 

A quantitative upper bound that immediately follows from ()3.4()j) is 



0< \(S Pt ,5q t )\ <s9t\\S e t £ \\ ns sup ||Sf|| H S7 I \V'(p e /,q e /)\ds 

se[o,t] J Jo 

<etm) 2 sup i f\V\p e ^q d f)\ds. 

ee[o,i] 1 Jo 



(3.41) 



0e[o,i] 

Here we have introduced 

Sf= sup ||Sf||HS, (3.42) 

se[o,t],6»e[o,i] 

whose asymptotic behaviour in the case of linearly unstable trajectories is controlled by 
the exponent 

T = sup X d£ . (3.43) 
ee[o,i] 
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Furthermore, under favourable circumstances the time average V of V in ()3.41|) is finite 
as t — > oo; then the asymptotic behaviour of the r.h.s. in (j3.41|) for large times is given by 

eV } te 2Tt . (3.44) 

This will, e.g., be the case if either the trajectory remains in a bounded set, or the derivative 
V is a bounded function on the respective energy shell. 



4 Inverted oscillators 

We want to discuss a simple and exactly solvable example that nevertheless possesses 
the typical features of exponentially unstable classical dynamics: a rf-dimensional inverted 
harmonic oscillator. In that case the Hamiltonian is quadratic in position and momentum 
and therefore the semiclassical propagation ()2.7|) is exact, i.e., the error term vanishes. 
To be specific, let 

H (p 1 q)= 1 -p 2 -^-q\ (4.1) 

and define 

H £ (p,q) = H (p,q-ea) , (4.2) 

so that, up to a constant, V(q) = to 2 a ■ q. This perturbation consists of a phase space 
translation of the unperturbed Hamiltonian and hence is of the same type as the one 
discussed in |20 a . The equations of motion generated by the unperturbed and by the 
perturbed Hamiltonian, respectively, can be solved explicitly, leading to 

p? = p cosh tot + quj sinh tot , 

_ (4.3) 
q° t = q cosh tot + p to' 1 sinh tot , 



and 



5p t = aeto sinh tot , 

/ . (4.4) 

5q t = ae(coshtot — 1) . 

From (|2.1()j) . ()4.3j) . and f)4.4j) one, moreover, obtains 

gs _ go _ f cosh^tl uosmhuot l\ ,^ ^ 

* ' v^^ 1 sinh uot 1 cosh tot 1 J 

This implies an accidental coincidence of the unperturbed and the perturbed Lyapunov 
exponents: X° = to = X £ . Furthermore, Ggo^ = Gg*[z] so that -Fdi sp (t) = 2~ d , reflecting 
the fact that the coinciding perturbed and unperturbed linearised dynamics lead to the 
same dispersions of the wave packets. Since 5X = 0, this finding is in accordance with the 
bounds ()3.23|) . With the help of the relation (|A.8|) and a change of variables the expression 
(13. 3 j) for the fidelity can now be brought into the form 

F(t) = JI Wmii^y) - S-\5p t ,5q t )) Wm^y) |^ . (4.6) 
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Therefore, in this example the quantum fidelity is crucially determined by the separation 
()4.4|) of the classical trajectories. 

For simplicity one can imagine the initial wave packet to be a Gaussian (|2.1J) with 
Z = it that is localised at the unstable fixed point (p, q) = (0, 0) of the unperturbed 
classical dynamics. Thus, q°) = (0, 0) for all t, so that the unperturbed time evolution 
()2.7j) fixes the centre of the wave packet and only forces it to disperse according to the 
action of the metaplectic operator related to ()4.5|) . The perturbed dynamics, however, 
pushes the centre away from the fixed point according to ()4.4|) . This happens with an 
exponential rate that follows from the asymptotic behaviour 

\(5p t ,5q t )\~ £ ^VvTTe xt , t^oc, (4.7) 

of the distance, which may be compared with the corresponding asymptotics 

£|a|A 2 te 2A * (4.8) 

of the upper bound (j3.41|) . see also ()3.44|) . 

For Gaussian states either (j3.12|) or ()4.6|) can be evaluated directly, yielding 

2 2 

F(t) =exp -^((l-coshAt) 2 -A 2 sinh 2 At) . (4.9) 

In this example the quantum fidelity therefore decays extremely fast, namely in a double 
exponential manner. We stress that neither have approximations been performed nor have 
any assumptions entered, and hence the result holds unconditionally. Clearly, this finding 
is at variance with the previous predictions of an exponential decay of fidelity. However, 
it obviously complies with the bounds ()3.29|) and ()3.34|) derived in Section EJ We remark 
that the double exponential decay is caused by both the separation ()4.4|) of the trajectories 
and the exponential instability of the linearised dynamics (|4.5jl . Each of these factors alone 
would lead to this effect. 



5 Conclusions 

Our approach to the quantum fidelity of localised wave packets led us to distinguish two 
effects that derive from two separate contributions to the semiclassical evolution prior to 
the Ehrenfest time. 

One effect is caused by the dispersion of the wave packets, which semiclassically origi- 
nates from the metaplectic representation of the linearised classical dynamics. Since gener- 
ically the unperturbed and the perturbed classical dynamics possess different linearisa- 
tions, the resulting non- coinciding dispersions cause an eventually exponential contribution 
Fdis P (t) to the decay of fidelity as described by ()3.22|) and ()3.23|) . 

A second effect is due to the separation of the perturbed classical trajectory away from 
the unperturbed one. Since the centres of the wave packets follow their associated classical 
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trajectories, this divergence forces the overlap of the unperturbed with the perturbed time 
evolution of the given initial state to decrease. We estimated this contribution, -F class (i), 
for Gaussian wave packets and identified an influence of the linear stabilities as well as of 
the separation of the trajectories. Since in general the latter cannot be well controlled, we 
were unable to determine a uniform expression for this factor. The bounds (|3.29jl . ()3.34j) . 
and (j3.41|) that we obtained allow for decays that are exponentially faster, or slower, than 
exponential. Of course, the further factor F disp (t) always ensures that the fidelity decays 
at least exponentially. 

In view of the previous predictions of an exponential fidelity decay a contribution that 
decreases in a double exponential manner might come as a surprise. In the example of 
inverted harmonic oscillators, however, we saw that such a behaviour is indeed possible. 
In that case this was caused by both the linear instability of the classical motion and by 
the exponentially growing separation of the trajectories. The latter effect is certainly not 
generic if one has chaotic systems with bounded energy shells in mind. Nevertheless, in our 
example the linear instability alone would cause a doubly exponentially decreasing factor. 
And this is in perfect agreement with the bound ()3.29|) that applies in the general case, 
even if the separation of trajectories does not exceed a given bound. 
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Appendices 

A Metaplectic representation 

In this appendix we collect some important facts about the symplectic group and the 
metaplectic representation. For further details see ITol ITTj. 

The symplectic group consists of the linear canonical transformations (p, q) h- > (p', q') 
with 



Ap + Bq , 
Cp + Dq . 



(A.l) 



The real 2d x 2d matrix 




(A.2) 



then fulfills S T JS = J, where 





(A.3) 
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is the symplectic unit. The symplectic group is generated by the matrices 



s -< = (o (AO ■ Sc = (5 3 ■ J ■ (A4) 

where A is an invertible matrix and C is symmetric. 

A quantisation of a linear canonical transformation (|A.1J) requires a unitary ray-re- 
presentation of the symplectic group. This can be obtained from the observation that 
the operators D(p,q) and Z)((5' T ) _1 (p, q)), see (|2.4|) . each provide a unitary irreducible 
representation of the Heisenberg group. According to the Stone- Von Neumann Theorem 
there hence exists a unitary operator £l(S) such that 

D({£f r )- 1 {p,q)) =iM(S)D(p,q)fi(S)- 1 . (A.5) 

Choosing 5 = S1S2 furthermore implies the multiplicative property 

ftSxSs) = e'x&M KSi) ■ (A.6) 

In fact, the phase factor can be chosen to be ±1. The metaplectic operators p.(S) deter- 
mine a double-valued representation of the symplectic group which is also known as the 
metaplectic representation. 

Up to a sign the metaplectic representation is fixed once the metaplectic operators for 
the generators (|A.4|) are given. Exploiting the relation (|A.5[) . one obtains 



fi{S A )ip{x) = Vdet A^j(A t x) , 

jl(S c )ij(x) = ±e^ x - Cx ij(x) , (A.7) 

/i(j)V>(p) = i d/2 V<p) , 

for them, where ip(p) denotes the momentum representation of ip. 

An explicit calculation based on the relation (|A.5|) finally reveals that the Wigner 
representation of a quantum state is covariant under linear canonical transformations, 

W\fi(S)if>](£,x) = W[il>]{S-\t,x)) . (A.8) 

Thus, if ip is localised at the point (p, q) in phase space, the transformed state fi{S)ip is 
concentrated at S(p,q). 



B Diagonal form of positive matrices 

It is well known that if two real and symmetric matrices commute, they can be simultane- 
ously diagonalised by an orthogonal transformation. Less appreciated is the possibility of 
converting non-commuting, but positive-definite, symmetric matrices into a diagonal form 
with a single transformation: 
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Let A and B be real, symmetric, and positive-definite n x n matrices. Then there 
exists a real, invertible (not necessarily orthogonal) matrix M such that M T AM = 1 and 
M T BM is diagonal, with the eigenvalues Aj(A~ 1 B) of the positive-definite matrix A~ X B 
on the diagonal. Moreover, 

n 

det{A + B) = Y[A J (A)[1 + A ] (A- 1 B)} . (B.l) 
i=i 

A proof of this statement is not difficult: Let O be an orthogonal matrix such that T AO = 
D is diagonal (and positive-definite). Define U = OD~ l l 2 , then U T AU = 1, and U T BU 
is symmetric and positive-definite. Furthermore, since U T BU = U^A-^BU, the matrices 
U T BU and A~ 1 B have identical spectra. Hence there exists an orthogonal matrix 0\ such 
that OiU T BUOi is diagonal, with the eigenvalues Aj(A~ 1 B) on the diagonal. Then define 
M = UOi to obtain the matrix M of the statement. 



C Matrix norms and singular values 

Real n x n matrices can be estimated in terms of various norms, for which there exists a 
number of inequalities that we want to review in this appendix. More details can, e.g., be 
found in fH\ . 

The operator norm is defined as 

1 1 ^4 1 1 op = sup \Ax\ , (C.l) 
l*l=i 

where \x\ = \fx? is the euclidean norm of a vector x 6 W 1 . The trace norm, however, is 

given by 

\\A\\ tv = ti VA^A . (C.2) 

Finally, we consider the Hilbert- Schmidt norm 

\\A\\ HS = VtxA^A , (C.3) 

All of these matrix norms possess the multiplicative property ||AB|| < ||A|| Moreover, 
they fulfill 

\\A\\o P < \\A\\us < \\A\\ tI , 
||^-B||tr < ||A||hs||-B||hs • 

In addition, for symplectic matrices S one obtains ||5' _1 || tr /Hs = || -S'Htr/HS- 

In general a real n x n matrix A cannot be diagonalised. However, A 7 A is symmetric 
and positive-definite and therefore possesses n non-negative eigenvalues. Their positive 
square- roots ^j{A) are the singular values of A, which we order as 

/U max (A) = fii(A) > fi 2 (A) >■■■> n n {A) > . (C.5) 
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Furthermore, Fan's inequality (see [21]) implies for the singular values of products that 



/i k (AB) < /w(#) fx k (A) . [ ' ' 

The singular values of symplectic matrices occur in pairs of mutually inverse numbers. 
Since in that case n must be even, we write n = 2d, and choose the following ordering: 

Hi(S)>---> fx d {S) > ^(S)- 1 >■■■> ^(S)- 1 . (C.7) 



In addition to the upper bound (|C.6[) . in the symplectic case one can also find a lower 
bound that is based on the fact that /imax(5'~ 1 ) = / i max(5')- Choose first A = SiS 2 and 
B = S^ 1 , and then A = S± 1 and B — S1S2 in (|(L6p . This results in 

max ( ^-4^rr , ^-tttt 1 < ^k{SiS 2 ) < min{/i fc (5'i)/i max (S'2),/ifc(S'2)/i max (5'i)} . 



(Si) 

(C.8) 

In section 13.21 we need to estimate a quadratic form Bv ■ ABv in terms of suitable norms 
of the positive-definite, symmetric matrix A and of the invertible matrix B. An upper 
bound follows from the definition (jC.lj) of the operator norm and a subsequent application 

of (EH, 



Bv ABv < ||A|| op ||5pw| 2 



< \\A\\ tT \\B\\ 2 m \v 



2 



(C.9) 



A lower bound can be gained from the fact that the quadratic form defined by A attains 
its minimum at the eigenvector corresponding to the lowest eigenvalue A m i n (A) > 0. Thus 



Bv ■ ABv > A min (A) \Bv\ 2 

= A min (A)v-B T Bv (CIO) 
> A min (A) /i n (B) 2 \v\ 2 , 

since by (1C.5J) fi n (B) is the lowest singular value of B. 
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